Unraveling the Nature of Charge Excitations in La2Cu04 with Momentum- Resolved 

Cu /\-edge Resonant Inelastic X-ray Scattering 
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Results of model calculations using exact diagonalization reveal the orbital character of states 
associated with different Raman loss peaks in Cu K-eAge resonant inelastic X-ray scattering (RIXS) 
from La2Cu04. The model includes electronic orbitals necessary to highlight non-local Zhang- Rice 
singlet, charge transfer and d-d excitations, as well as states with apical oxygen 2pz character. The 
dispersion of these excitations is discussed with prospects for resonant final state wave-function 
mapping. A good agreement with experiments emphasizes the substantial multi-orbital character 
of RIXS profiles in the energy transfer range 1-6 eV. 



PACS numbers: 78.70.Ck, 78.20.Bh, 74.25. Jb, 74.72.Cj 

Understanding the nature of charge excitations in cor- 
related systems remains a challenge even after more than 
two decades of intense study. Recent experiments in the 
cupratc high-Tc superconductors highlight the potential 
importance of additional orbital degrees of freedom in- 
volving the apex oxygens. [1] Moreover, a correlation be- 
tween the apex oxygen site and Tc has been made in 
different contexts. [2] However, to date no experiment 
has been able to directly test the apical character of the 
low energy states in the cuprates. Nonetheless, these re- 
sults suggest a re-examination of the canonical view that 
the physics in cuprates is solely governed by electrons 
confined to the in-plane copper-oxygen orbitals. 

One paradigmatic material is La2Cu04, on which 
many spectroscopic techniques have been used to elu- 
cidate various aspects of the underlying physics. [3] 
Among these techniques, resonant inelastic x-ray scatter- 
ing (RIXS) [4] probes charge excitations in a momentum- 
resolved way. A net momentum AQ is transferred during 
the RIXS process, and the associated photon wavelength 
allows access to the entire Brillouin zone (BZ). By tuning 
the x-ray energy to proper absorption edges, RIXS pro- 
vides material-specific information of various many-body 
excitations, e.g. magnon, d-d, and charge transfer (CT) 
excitations, in strongly correlated materials. [4-9] 

Cu K-edge (ls-4p core- level excitation) RIXS measure- 
ments in the cuprates reveal a number of spectral peaks 
in the energy-loss range 1-6 eV, [5-9] with a most promi- 
nent peak ~ 4-6 eV and another less prominent peak '^2 



eV at AQ = (0,0) (the T point or zone center). The 
former has been attributed to a local CT of a hole from 
Cu to its neighboring oxygen ligands, and the latter has 
been attributed to excitation from the lower Hubbard 
band across the band gap, leading to a CT excitation ac- 
cording to the Zaancn-Sawatzky- Allen scheme. [10] These 
and other excitations have been studied carefully in effec- 
tive low energy Hamiltonians. [11] However, these down- 
folded approaches fail to fully account for the impact of 
many intertwined orbitals on RIXS profiles. 

In this letter, a "two-step" exact diagonalization (ED) 
algorithm facilitates the characterization of the Cu K- 
cdgc RIXS spectra in La2Cu04 in the energy range 1-6 
cV with various AQ. These excitations is shown to be 
rich in orbital character and not described by the dis- 
persion of a single Mott gap excitation. The d-d, apical 
oxygen, and CT excitations are all found in the studied 
energy loss range, having intensities that vary with AQ. 
The method consists of (1) a matrix diagonalization of 
the model Hamiltonian to obtain the ground state and 
final excited states, and (2) a series of diagonalizations to 
obtain intermediate states with the core-hole interaction. 
We use the parallel version of the ARPACK libraries, 
based on iterative Arnoldi methods, [12] to calculate the 
energy-eigenvalues and orthonormalized eigenfunctions. 

We focused on quasi-2D copper-oxygen clusters with 
periodic boundary conditions. The CU4O2; cluster allows 
investigations of three distinct AQ-points: (0,0), (7r,0) 
and (7r,7r), which is enough to demonstrate the qualita- 
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FIG. 1: (Color online) Theoretical RIXS spectra at the F 
point for (a) planar CU4O8, and (b) planar CugOie cluster 
(with the elastic resonance removed). The red dashed rectan- 
gle in (a) indicates the energy range of the RIXS plot in (b) 
(intensity xlO^). The spectra are broadened with a 0.1 eV 
Lorentzian in energy loss and a core hole lifetime of 1.0 eV. 



tive differences in tfie momentum-dependent RIXS spec- 
tra. We use a multi-orbital Hubbard Hamiltonian includ- 
ing Cu 3d Eg and bonding O 2p orbitals, with the Hund's 
coupling and direct exchange interaction to account for 
multiplct splitting. [13, 14] Explicit photon-polarization 
dependence is introduced by convolving the RIXS spec- 
trum with a Cu 4p DOS obtained from density functional 
theory (DFT) calculation. [15] 

Results at the F point are displayed in Figs. 1 and 
2, where we show progressively more spectral complexity 
by including more orbitals. Retaining only Cu Sdx^~y^ 
and its hybridized planar O 2px,y orbitals [CU4O8] yields 
Fig. 1(a). Besides the elastic feature, there is only a 
single prominent peak at energy transfer ^ 4.8 eV. This 
peak stems from a local CT excitation of a Cu cP hole 
onto its surrounding ligand oxygens with d^^-y^ symme- 
try, creating a local molecular orbital d^^L. This energy 
loss is directly controlled by planar charge transfer energy 
A = — Erf, and is further increased by a level repulsion 
with the ground state via p-d hybridization. The inclu- 
sion of an excitonic interaction between the Cu Zd elec- 
tron and O 2p hole would tend to decrease the excitation 
energy. However, for the ligand state of dxi_y2 symmetry 
this effective electron-hole interaction is small. [16, 17] 

Wc note that while the energy gap revealed by the cal- 
culated hole addition/removal spectra is '^1.7 eV, [18] 
there is no RIXS resonance at this energy at the F point 
associated with a Zhang-Rice singlet (ZRS) CT excita- 
tion. Similar calculations have been performed for a pla- 
nar CugOie cluster up to energy loss 2.2 eV. As shown 
in Fig. 1(b), a weak 0.56 eV bi-magnon peak is found, 
but a final state resonance '^2 eV remains absent. 

As shown previously, [19] a non-local ZRS CT exci- 
tation is well-defined away from the F point (especially 
along the magnetic BZ boundary). Although the 2 eV 
ZRS feature is absent in the hole addition spectra at the 
F point measured at half-filling by angle-resolved pho- 
toemission due to complete cancellation between 2px^y- 
3dx2-y2 hybridizations, [20, 21] its absence in RIXS can 
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FIG. 2: (Color online) Theoretical RIXS spectra at the F 
point for (a) octahedral CU4O16 with apical O, (b) planar 
CuJOg (with Cu 3d322_r2), and (c) octahedral CuJOie (with 
both Cu 3d^^2_^2 and apical O). (d)-(f) Schematic of the 
non-local inter-site d-d excitation pathway for the results in 
(b) and (c): (d) ground state configuration, (e) intermedi- 
ate ligand occupation on both plaquettes, and (f) final state 
configuration highlighting the inter-site d-d excitation. 



be more subtle. Unlike ARPES, in A'-edge RIXS the hole 
number is conserved, and a d^L is created together with 
a via charge transfer. Certain configurations, where 
the momentum of a ZRS and a d^*^ hole sum to zero, may 
admit a final state resonance associated with ZRS at the 
F point. Therefore, the absence of a 2 eV feature may be 
an effect of confinement on the d^L and d^" and thus an 
artifact of the small cluster size. To provide a definitive 
answer regarding the existence of a 2eV feature in RIXS, 
it is essential to go beyond the CugOie calculation. 

We show results at the F point that include apical oxy- 
gen 2pz, copper 3d322_^2, and both orbitals in Fig. 2 (a)- 
(c), respectively. Apical oxygens [CU4O16 in Fig. 2(a)] 
yield an additional peak compared to Fig. 1 at an en- 
ergy loss ~ 3.5 eV controlled by the apical CT energy, 
attributed to a nonlocal symmetry allowed CT excita- 
tion of the Cu d^ hole onto inter-plaquette apex oxygens. 
Although the CT energies A are similar for planar and 
apical oxygen orbitals, this apical excitation lies below 
the d^°L due to a small level repulsion with the ground 
state, stemming from a lack of 2pz-3d^2_y2 hybridization. 

Inclusion of the Cu id'^y2_^2 orbital (Figs. 2(b) 
[CulOg] and 2(c) [Cu|0i6]) produces peaks at lower en- 
ergies. These low energy peaks, present only with the 
d3z2_^2 orbital, are of d-d character. An intra-site d-d 
excitation, where a Cu 'id^z^_j.2 hole is left behind after 
Ap-\s recombination on the same Cu atom, is symmetry- 
forbidden. [9, 15] However, non-local inter-site d-d ex- 
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FIG. 3: (Color online) Theoretical RIXS spectra for CU4O8 
at AQ = (a) (0,0), (b) (tt.O), and (c) (tt.tt). (d), (e), and 
(f) represent the same momentum points for CU4O16. The 
dispersions of local molecular orbital and inter-plaquette api- 
cal CT excitations are indicated respectively by the blue and 
green dashed lines. 



citation, where a Cu 8^322 hole is excited from its 
neighboring Cu 3d^2_y2 orbitals, is symmetry-allowed 
[sec Fig. 2(d)-(f)]. The resonant features at loss ener- 
gies ~ 1.5 and 2.0 cV are characterized by single and 
double nonlocal d-d excitations. While the bare crystal- 
field-splitting (CFS) between the two eg orbitals in the 
calculation is only 0.2 eV, the de-localization between the 
Cu 3d^2_y2 hole and the ligand (^ 1.39 eV) [23] causes 
the peak to appear at higher energies, essentially the 
same effect as the covalent contribution to the ligand field 
splitting. [22] More peaks emerge around the d^'^L molec- 
ular orbital excitation which contain different mixings of 
Cu 3d^2_y2 ,3d^z2_,^2 and O 2px,y,2pz orbitals. More- 
over, Figs. 2(b) [Cu^Os, relevant for Sr2Cu02Cl2 and 
Nd2Cu04] and 2(c) indicate that the presence of apex 
oxygens shifts the RIXS peaks to lower energy, due to 
kinetic energy stabilization of states with Sd^z^-r^ char- 
acter. Taken as a whole, these results demonstrate that 
the orbital character of charge excitations accessible via 
RIXS becomes increasingly more complex and material 
dependent at higher energy transfer. 

In Fig. 3 we describe the dispersion of various RIXS ex- 
citations. First, for CU4O8 [Figs. 3 (a)-(c)] the prominent 
molecular orbital excitation energy softens from (0, 0) to 
(tt, 0), and increases again at {tt,tt) [follow the dashed 
line in Fig. 3(a)], yet remains lower than its energy at 
the r point, indicating an indirect energy gap. This be- 
havior agrees with experimental observations. [5] At fi- 
nite momentum transfer AQ — (tt, 0), d^'^L excitation 
propagates more easily because of its associated ligand 
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FIG. 4: (Color online) (a) Experimental [6] and (b) theo- 
retical RIXS spectra for different AQ at incident photon en- 
ergy 8992.5 eV. The elastic part for AQ = (0, 0) is removed. 
The spectra in (b) are broadened with an energy-dependent 
Lorentzian increasing linearly with energy loss from 0.1 to 0.4 
eV, 0.2 to 0.5 eV, and 0.3 to 0.6 eV, for AQ = (0,0), {tt,0), 
and {n,n), respectively, (c) Experimental and (d) theoretical 
RIXS spectra convolved with 4pz DOS at AQ = (7r,7r). 



symmetry. As shown previously, [19] a non-local ZRS CT 
excitation is well-defined along the magnetic BZ bound- 
ary and does not suffer from the potential impact of the 
small cluster size. This is confirmed in Figs. 3 (b) and 
(e), where a smaller non-local ZRS CT peak emerges at 
2.8 eV at AQ = (7r,0), --800 meV higher in energy 
than that of the anticipated ZRS at the F point. 

The Cu|0i6 dispersion profile [Fig. 3(d)-(f)] can be 
understood with help from the CU4O8 analysis. The 
dominant peak at AQ = (tt, 0) at 3.5 eV results from 
the local molecular orbital excitation. On the other hand, 
the ~ 3.5 eV inter-plaquette apical CT excitation at the F 
point has little mixing with the ground state, and there- 
fore only weakly disperses. The RIXS data, focusing on 
the dispersion of the excited states and their characters, 
thereby provide a vehicle to probe the apical character 
of the many-body states. The cluster calculations give 
only a few percent apical weight in the ground state, [24] 
implying that phenomena associated with apical ground 
state character could be difficult to detect. 

Figure 4(a) displays cuts of experimental RIXS spectra 
on La2Cu04 at three different AQ points. [6] The RIXS 
profile broadens into a single peak as AQ increases. As 
higher energy states lie in a high density continua, we 
use energy- dependent broadening for different AQ. A 
generally good agreement with the experiments at finite 
AQ results from this procedure, particularly for energies 
above 3 eV. However, the calculation does not obtain 
the broad intensity between energy loss 2 3 eV at the 
F point, and seems to overestimate the d-d excitation 
intensity. Agreement may be improved at lower ener- 
gies with larger cluster studies or smaller Hilbert space 
calculations with downfolded Hamiltonians, which would 
emphasize the ZRS CT with a likely suppression of d-d 
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spectral weight. We also remark that a nonlocal ZRS 
CT excitation with momentum AQ may appear at the 
r point with another excitation, such as a phonon or a 
magnon, carrying momentum — AQ. 

Finally we demonstrate that incorporating 4^^ DOS 
obtained from DFT produces an adequate position for 
the K-edge resonance and the main energy-loss feature. 
RIXS data for La2Cu04 from Ref. 8 are replotted in Fig. 
4(c). The spectra were collected at room-temperature in 
vertical scattering geometry, with e || c (AQ=(2.5,0.5,0)) 
and overall energy resolution set to ^ 0.3 cV. Experimen- 
tal RIXS spectra have broad features, with a dominant 
peak at incident photon energy ^ 8992 eV and energy 
loss ^ 4 eV, agreeing with the calculation [Fig. 4(d)]. 

In summary, this paper describes a nontrivial momen- 
tum dependence of the RIXS process and the resulting 
multi-orbital character of the RIXS profile for La2Cu04. 
We demonstrate how the multi-orbital character of ex- 
cited states may be mapped throughout the BZ, and give 
a reason why a dispersion of certain peaks may be asso- 
ciated with a change of character. Moreover, we show 
that subdominant character of the ground state may be 
revealed as apparent dispersion of a local excitation. 

We also demonstrate how the presence of the apical 
oxygens can affect the ground state and the excitation 
spectrum. From the spectral weight and dispersion, one 
can confirm whether there is substantial apical character 
of the ground state. As the apical oxygen energy is one of 
the few parameters that changes across the high Tc com- 
pounds, K-edge RIXS in the energy range 1-6 eV may 
offer a diagnostic tool to examine material dependence, 
and provide hints to the material-specific Tc. [25] 

For the undopcd cluster, the apical oxygen character, 
including that associated with a possible time-reversal 
symmetry breaking (TRSB) phase involving circulating 
current loops through the apical oxygens, has very small 
weight in the ground state, with the majority appear- 
ing in higher energy excited states. This is consistent 
with variational Quantum Monte Carlo and cluster t-J 
results which found evidence for such a phase only for 
larger doped clusters. [26]. RIXS in principle can shed 
light onto a putative TRSB phase by performing circu- 
larly polarized dichroism experiments at the correspond- 
ing ordering wave- vector. An extension of the current 
calculations to address this interesting point for doped 
multiband clusters is an area of future study. 
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